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SUMMARY 


Theoretical  expressions  are  derived  for  the  aerodynanuc  forces 
on  a  helicopter  rotor,  and  for  their  contributions  to  various  stability 
derivatives.  The  resulting  expressions  are  tabulated,  containing  a 
number  of  functions  of  advance  ratio  which  are  presented  graphically. 
Offset  hinges  and  pitch-flap  coupling  are  provided  for  in  the  formula¬ 
tion.  A  novel  method  for  trim  calculation  is  presented. 

The  theory  is  based  on  the  original  results  of  Bailey  (Ref.  13). 
It  involves  some  changes  in  notation,  axes,  assumptions,  and  the 
structure  of  formulas  to  widen  their  scope  and  improve  their  con¬ 
venience  for  stability  estimates.  Discussion  of  the  limiting  assump¬ 
tions  is  offered. 
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1 .  Introduction 

The  material  of  this  report  was  originally  composed  by  the  authors  for 
a  handbook  on  helicopter  stability  and  control,  sponsored  by  the  U.S.  Navy. 
Although  the  handbook  was  subsequently  published  (Reference  16),  its  dis¬ 
tribution  and  use  have  not  been  extensive,  and  it  was  felt  that  the  separate 
publication  of  the  chapter  on  rotor  aerodynamics  would  be  valuable  as  a  basis 
for  VTOL  as  well  as  helicopter  analyses. 

In  this  report,  certain  minor  errors  of  the  original  have  been  corrected 
and  the  structure  of  some  of  the  formulas  has  been  simplified.  The  brief 
treatment  of  fuselage,  tail  and  tail  rou>r  stability  contributions  seems  now 
too  incomplete  to  be  very  useful,  and  so  it  has  been  deleted.  The  considerable 
effort  in  the  original  work  to  present  rotor  aerodynamics  in  terms  as  general  as 
possible  has  been  preserved,  so  that  the  equations  might  be  useful  for  heli¬ 
copters  of  different  configurations,  with  rotors  of  at  least  some  different  im¬ 
portant  types . 

There  are  Ftill  many  areas  of  incomplete  understanding  in  rotor  aero¬ 
dynamics.  Much  of  the  work  presented  here  rests  on  still  doubtful  assump¬ 
tions,  and  the  results  are  still  largely  uncnecked  by  experiment.  The  working 
assumptions  and  their  probable  limitations  are  therefore  to  be  carefully  noted. 
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2 .  Notation 

2.1.  forces,  Moments,  and  Axes 

In  this  report,  rotor  forces  and  moments  acting  on  the  aircraft  are 
measured  along  and  about  three  mutually  perpendicular  axes,  of  which  the  H 
axis  Is  parallel  to  the  rotor  drive  shaft,  and  the  X  lies  in  the  fuselage  sym¬ 
metry  plane.  For  these  rotor  forces  and  moments,  the  origin  is  the  rotor  hub, 
and  subscript  (  identifies  resolution  along  these  shaft,  axes.  Thus: 

Ts  Rotor  thrust,  along  shaft,  positive  upwards,  lbs. 

Hg  Rotor  "horizontal"  force,  J-  shaft  and  in  fuselage  plane 

of  symmetry,  positive  to  rear,  lbs. 


Rotor  "side"  .force,  J.  plane  of  Tg  and  Hg  positive  in 

0 

direction  of  iff  =90  (to  right,  for  counterclockwise 
rotation),  lbs. 


I  • 

Rotor  force  vector  tilt*  radians:  ** 

Rotor  hub  pitching  moment  about  axis  of 
nose  up,  ft-lbs. 


Hr 

Ys  ,  positive 


Lg  Rotor  hub  rolling  moment  about  axis  of  Hg  ,  positive 
right  side  down  for  ccw  rotation,  ft-lbs. 


Q  Rotor  torque  about  shaft  axis,  positive  decelerating,  ft-lbs. 

Rotor  "wind"  axes  are  also  used  to  find  formulae  for  the  "shaft"  forces  and 
moments  above.  The  “wind"  axes  are  mutually  perpendicular,  with  origin  at 
the  rotor  hub.  The  2y/  axis  is  along  the  axis  of  no  feathering  and  the 
axis  lies  in  the  plane  of  Zw  and  the  relative  wind  at  the  rotor  hub.  Thus: 


0. 


Tw  Rotor  thrust,  along  ANF  ,  positive  up,  lbs. 

Hw  Rotor  "horizontal"  force,  JL  ANF  ,  in  plane  of  Tw 

and  relative  wind,  positive  rearwards,  lbs. 

Yw  Rotor  “side"  force,  _L  plane  of  Tw  and  , 

positive  in  direction  of  ^  *90° (to  right  for  CCW 
rotation),  lbs. 

Mw  Rotor  hub  pitching  moment,  about  axis  of  Yw  ' 

positive  nose  up,  ft-lbs. 


t  Rotor  hub  rolling  moment,  about  axis  of  Hw  ,  positive 

right  side  down  for  CCW  rotation,  ft-lbs. 

All  forces  are  rendered  non-dimensional  by  dividing  by  />IttF^XQK)Z. 


Thus: 


CT  S  r  2 
*  pirRZ  (S  R)z 

2  2 

All  moments  are  rendered  non-dimensional  by  )( £l  R)  R/ 


/>irR*  (OR)2 


Thus: 
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2 . 2  Aircraft  and  Rotor  Physical  Characteristics 


w  Total  aircraft  weight,  lbs. 


b  No.  of  blades. 


R  Rotor  radius,  ft. 


cr  Solidity  ratio,  CT  = 


C  Equivalent  blade  chord,  ft. 


Lock's  blade  inertia  constant. 


/»ocR 


Blade  moment  of  inertia  about  the  flapping  hinge,  slug-ft 


Mass  moment  of  blade  about  flapping  hinge,  slug-ft. 


b  Flapping  hinge  eccentricity,  ft. 


Coupling  between  blade  pitch  and  blade  flapping,  due  to 

x  ~~  a/Sty) 

(positive  when  upward  flapping  reduces  blade  pitch) 


I 
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*f|  i^2r*4  Functi°ns  of  r  and  y  ,  see  Figures  8  to  19. 

f ! ,  fg%...  Functions  of  r  and  y  ,  ana  related  to  above  f  j  , 
fg  ,  etc.  by 

f'r.  •  +„(T)-t#(0) 


2-3.  Aircraft  Flight  Condition 

V  Velocity  along  the  flight  path*  ft.  per  sec. 

€  Climb  angle,  positive  in  climb,  radians. 


Angle-of-attack,  between  the  remote  wind  vector  and  a 
plane  JL  rotor  shaft,  positive  nose  up,  radians. 


Sideslip  angle,  between  remote  wind  and  plane  of  sym¬ 
metry,  positive  when  wind  is  from  the  right,  radians. 


Aircraft  ra<"5  of  yaw  about  2$  axis,  positive  nose  right, 
radians/sec.  (Aircraft  yaw  angle  does  not  occur  in  the 
equations.  Symbol  \f/  is  used  for  blade  azimuth 
angle* 


ps  Aircraft  roll  rate  about  Xg  axis,  positive  in  right  roll, 
radlans/sec. 

<S  Aircraft  pitch  rate  about  Ys  axis*  positive  nose  up, 
radians/sec. 


i 

i  ' 
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U§  Ratio  of  Xg  velocity  component  to  tip  speed 

_  V  cos  as 
“s=  “SR 


Wg  Ratio  of  Zg  velocity  component  to  tip  speed 

_  _  V  gin  a8 
ws  -  £  R 


2.4.  Rotor  Operating  Conditions  and  Aerodynamics 


a 


Advance  ratio, 


V 

flR 


Rotor  angular  velocity,  rad/sec.  (All  signs  in  force  and 
moment  equations  in  this  report  are  for  CCW  rotation 
of  main  rotor.) 


Ctyy  Rotor  angle  of  attack,  between  remote  win«  at  rotor  hub 
and  plane  JL  ANF,  positive  nose  up,  radians 
unless  noted  —  see  equation  18.  In  X  partial 
derivatives,  sub  {  )w  is  omitted. 

Apparent  sideslip  angle  at  rotor  hub,  positive  with  wind 
from  right,  radians  —  see  general  case  equation  11. 

a  Slope  of  blade  lift  curve,  usually  assumed  0  *  5.73/rad. 

8  Blade  mean  profile  drag  coefficient.,  often  assumed  8s.OI25 
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Xw  Total  inflow  ratio,  with  Inflow  velocity  parallel  to  ANF  , 
using  average  (over  disc)  induced  velocity  —  see 
equation  7.  In  formulas,  sub  (  is  omitted. 


X 


* 


Inflow  parameter,  ratio 


X.  Measure  of  fore-and-aft  variation  of  induced  velocity  -- 
1  j 

ratio  of  -g*  total  variation  (considered  linear)  to  tip 
speed. 


Average  inflow  ratio  for  hovering  at  particular  thrust 
coefficient 


xoE 


2 


a,.a‘2  Parameters  of  the  rotor  thrust  equations  —  see  equations  71# 

T2 ,  73  . 

U*  Li"  LI*  LJ*  Li*  U1  u' 

M0»M0  ,f12  ,m3  ,ri4 

Parameters  of  the  rotor  H  forcr  equations  —  see  equa¬ 
tions  74,  75,  76,  77,  78,79. 

y  i  y  ••  y  •  y  I  yl  yl  yl 

Y0,Y0  ,YI  » '2  ,y3  '  y4  *y5 

Parameters  of  the  rotor  Y  force  equations  —  see  equa¬ 
tions  80 ,  81  ,  82  ,  83  ,  84  ,  85 


2.5.  Rotor  Blade  Pitch 
General  equation  is 

5U)  "*0~  Alscos  *  “Bls8ln  * 


Blade  azimuth  angle,  measured  in  direction  of  rotation, 
from  downwind. 

General  blade  pitch  angle  as  function  of  azimuth,  radians. 

General  collective  pitch ,  function  of  control  and  coning  — 
see  equation  15. 

Collective  pitch  due  to  control,  or  8q  for  PM*  0.  With 
no  pitch-flap  coupling,  e0  *ec  . 

Longitudinal  tilt  between  shaft  and  ANF  ,  positive  ANF 
forward  of  shaft,,  radians.  In  general,  function  of 
control  and  cyclic  flapping  —  see  equation  16. 

Contribution  of  control  to  Sj^  ,  or  for  PM  *  o . 
May  be  associated  with  longitudinal  control  and  swash- 
plate  tilt.  Positive  stick  forward,  radians. 

Lateral  tilt  between  shaft  and  ANF*  positive  ANF  tilted 
toward  Vf  *90*  *  radians.  In  general,  function  of 
control  and  cyclic  flapping  —  see  equation  17. 

Contribution  of  control  to  Ajs  '  or  A|g  for  PW*  o. 
May  be  associated  with  lateral  control  and  swashplate 
tilt.  Positive  stick  right  for  ccw  rotation,  radians. 
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2.6.  Rotor  Blade  Flapping 
General  equation  is 

PW  s  a0  “  a,  cos  f  -  b,  sin  $ 

(subscripts  {  )j  or  (  used  to  differentiate  between  shaft  and  AN  Fas 
reference  axis.) 

PW  Blade  flap  angle  as  a  function  of  azimuth,  positive  up, 
radians . 

Qq  Coning  angle,  radians. 

0  j  Longitudinal  cyclic  flapping  or  tip  path  plane  tilt,  positive 

for  tip  plane  tilted  back;,  radians. 

b|  Lateral  cyclic  flapping  or  tip  path  plane  tilt,  positive  for 

tip  plane  tilted  toward  ^  s  90° (right  for  CCW  rota¬ 
tion),  radians. 
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3.  The  Working  Assumptions 

No  completely  rational  analysis  of  helicopter  aerodynamics  has  ever 
been  stated.  The  flows  involved,  especially  in  forward  flight,  are  extremely 
complex;  the  aerodynamics  of  the  rotor  blades  are  unsteady  and  non-linear; 
the  number  of  degrees  of  freedom  is  large — to  mention  but  a  few  of  the  com¬ 
plications.  Consequently,  any  reasonable  analysis  must  rest  on  a  multitude 
of  simplifying  assumptions,  of  which  there  is  a  great  variety  in  the  literature. 
The  validity  of  the  different  simplifications  varies  greatly  between  different 
analyses  according  to  the  purpose  for  which  they  will  be  used.  Fortunately, 
for  the  purposes  of  predicting  Handling  Qualities,  Stability,  and  Control, 
many  of  the  simpler  assumptions  appear  to  be  permissible. 

3.1.  Rotor  Induced  Velocity 

The  induced  velocity  is  assumed  constant  over  the  rotor  disc,  and 
given  in  value  by  simple  momentum  theory.  This  simplification  has  been 
almost  universally  adopted  for  Stability  and  Control  purposes  on  the  basis 
that  exact  rotor  forces  are  unnecessary,  but  rather  derivatives  of  change  in 
rotor  forces  with  perturbations  from  trim  are  important.  These  derivatives  are 
not,  for  the  most  part,  sensitive  to  the  details  of  the  induced  flow.  The 
major  defects  of  this  assumption  can  be  expected  to  be  associated  with  the 
real  fore-and-aft  variation  of  downwash  in  forward  flight,  which,  with 
articulated  rotors,  primarily  affects  lateral  blade  flapping.  This  in  turn 
may  affect  lateral  stability  and  trim,  and  provide  some  coupling  of  longi¬ 
tudinal  and  lateral  motions — the  latter  usually  neglected  am-way.  In  the 
material  to  follow,  even  the  fore-and-aft  downwash  variation  is  ignored, 
except  as  noted  in  connection  with  lateral  flapping. 

The  steady-state  momentum  value  of  the  induced  velocity  is  assumed 
to  apply  continuously,  changing  with  thrust  and  the  aerodynamic  variables. 
This  quasi-steady  assumption  is  probably  reasonable  for  most  aircraft  motions 
which  are  slow  enough  that  the  flow  response  time  could  be  considered 
negligible . 
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It  should  be  noted  that  in  vertical  flight,  the  constant  (over  the  disc) 
induced  velocity  assumption  is  inadequate  for  accurate  performance  predic¬ 
tion,  or  for  autorotation  analysis.  It  is  similarly  unacceptable  for  vibration, 
blade  motion,  or  blade  loading  analyses.  It  is  even  too  crude  for  Stability 
and  Control  studies  of  rigid  rotors,  where  the  longitudinal  downwash  variation 
in  forward  flight  produces  considerable  pitching  moments. 

Extended  discussion  of  the  alternatives  to  these  assumptions  about 
induced  velocity  is  beyond  the  scope  of  this  work.  A  comprehensive  survey 
may  be  found  in  Reference  1.  References  2,  3,  and  4  may  also  be  helpful. 
Suffice  it  to  say  that  any  refinement  in  this  matter  can  be  made  only  at  the 
expense  of  considerable  increase  in  the  complexity  of  the  analysis. 

3.2.  Rotor  Blade  Characteristics 

i.  Blades  are  rectangular  and  untwisted.  Most  present  day  helicopter 
rotors  use  blades  of  only  moderate  twist  and  taper,  because  of  manufacturing 
difficulties.  For  stability  purposes,  these  may  be  replaced  by  blades  of 
constant  equivalent  pitch  and  chord,  taken  as  the  real  values  at  3/4  radius. 
Appreciable  errors  will  be  avoided  if  the  trim  variables  (a ngle-of -attack, 
collective  pitch,  etc.),  which  occur  in  the  stability  derivatives,  are  cal¬ 
culated  on  the  basis  of  the  same  assumption.  For  performance,  autorotation, 
or  blade  analyses,  the  assumption  is  completely  unacceptable. 

ii.  Tip-loss  factor  equals  unity.  The  very  considerable  saving  in  com¬ 
plexity  of  equations  resulting  from  this  simplification  is  its  justification. 
Allowance  fo;  tip-loss  should  be  made  by  using  an  equivalent  rotor  radius 
equal  to  the  usual  tip-loss  factor  times  the  real  radius.  As  in  (i)  above,  the 
same  assumption  should  be  used  for  the  calculation  of  the  trim  variables. 

iii.  Blades  are  infinitely  stiff  in  torsion  and  bending.  This  assumption 
is  valid  for  blades  designed  for  coincidence  of  center  of  gravity,  elastic 
axis,  and  aerodynamic  center.  Any  deviation  from  this  condition  will  have 
important  effects  on  stability  and  control  characteristics.  The  more 
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significant  effects  can  usually  be  treated  by  superposition  of  suitable  incre¬ 
ments  on  the  basic  force  and  moment  derivatives  given  here.  For  so-called 
"rigid,"  or  non-articulated  blades,  the  assumption  may  introduce  serious 
errors,  even  with  balanced  blades.  Needless  to  say,  the  assumption  must 
be  discarded  for  blade  stability  and  strength  calculations. 

iv.  Blade  lift  is  zero  in  Reverse  Flow  Region:-  There  are  three  conven¬ 
tional  methods  of  treating  the  reverse  flow  region,  that  is,  the  area  on  the 
retreating  side  of  the  iotor  where  the  resultant  velocity  flows  from  the  trail¬ 
ing  edge  to  the  leading  edge.  The  first  method  is  to  ignore  the  change  of 
flow  direction  in  this  region  entirely.  Since  the  local  lift  of  the  blade  is 
dependent  upon  (u  +|^*sln  iff  ,  the  sign  change  of(U+ ^Sin  iff)  does  not 
appear,  and  so  thin  treatment  implies  lift  in  the  reverse  flow  region  opposite 
to  the  actual  direction.  The  second  method  is  to  revise  the  airload  equations 
such  that- the  change  of  lift  direction  is  taken  into  account.  Both  of  these 
methods  assume  that  the  airfoil  has  the  same  properties  in  the  reverse  direc¬ 
tion  as  in  the  forward  direction.  The  third  method  is  to  assume  that  the  lift 
of  the  blade  is  zero  in  this  region.  For  values  of  —  2  U  2s  I  ,  the  values 
of  horizontal  force  and  torque  predicted  by  each  of  these  methods  do  not 
differ  significantly.  However,  for-|0<  significant  differences 

appear.  Since  the  relative  velocities  over  the  blade  are  small,  and  at  an 
advance  ratio  of  .35  this  region  is  about  3%  of  the  total  disc  area,  it  does 
not  seem  physically  reasonable  that  the  actual  forces  should  vary  signifi¬ 
cantly  from  the  zero  lift  assumption.  Therefore,  the  assumption  of  zero 
lift  in  this  region  has  been  adopted  for  computation  of  the  horizontal  force 
and  torque.  Lack  of  experimr-  ntal  evidence  prevents  a  more  exact  treatment 
of  this  region. 

The  effects  of  the  reverse  flow  region  on  thrust,  sideforce,  and 
flapping  motions  are  believed  to  be  small,  no  matter  which  treatment  is 
used,  and  therefore  the  reverse  flow  region  is  ignored  entiieiy  in  these 
quantities . 
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v.  Blade  stall  and  compressibility  are  disregarded.  The  effects  of  these 
assumptions  on  Stability  and  Control  characteristics  are  minor  compared  to 
their  consequences  for  performance,  vibration,  or  blade  strength  calculations. 
In  most  cases,  flight  envelopes  will  be  limited  by  the  latter  considerations  to 
conditions  where  the  effects  on  Stability  and  Control  are  small.  The  subject 
is  one  demanding  further  treatment,  for  results  thus  far  are  rather  qualitative. 
References  5,  6.  and  7  are  recommended  for  detailed  study. 

3.3.  Rotor  Blade  Motions 

i.  Harmonics  of  blade  motion  above  the  first  are  negligible.  For  pur¬ 
poses  of  predicting  overall  aircraft  motions,  this  simplification  is  quite  gene¬ 
rally  acceptable.  The  forces  and  moments  created  by  higher  harmonic  blade 
motions  are  of  high  frequency  and  alternating,  with  zero  "average"  values. 

The  aircraft  has  sufficient  inertia  that  it  cannot  respond  significantly  to  these 
excitations,  and  their  importance  is  confined  to  the  areas  of  vibration  and 
blade  strength  considerations.  In  some  cases,  due  to  non-linear  effects, 
higher  harmonic  flapping  can  influence  coning  and  first-harmonic  motions 
(and  hence  also  the  Stability  and  Control)  but  these  interaction  effects  are 
usually  small  and  safely  neglected  for  the  present  purposes. 

ii.  Transient  blade  motions  are  replaced  by  quasi  static  blade  motions. 
This  assumption  usually  associated  with  Dr.  Hohenemser  (Reference  8),  re¬ 
fers  to  the  manner  in  which  the  blades  follow  the  motion  of  the  fuselage.  The 
essential  concept  is  that  the  natural  frequency  and  damping  of  the  blade 
motion  with  respect  to  the  axis  of  no  feathering  are  so  high  that  the  blade 
can  be  considered  as  responding  instantaneously  to  fuselage  motions.  The 
blade  flapping  angle  at  any  instant  can  therefore  be  calculated  by  considera- 
tion  of  the  rotor  operating  conditions  at  that  instant,  where  pitch  rate  and 
roll  rate  must  be  included  in  the  variables  describing  the  condition.  The 
gyroscopic  forces  acting  on  the  blades  due  to  pitching  or  rolling  create 
steady  flapping  in  a  direction  to  resist  the  roll  or  pitch.  This  effect  results 
in  damping  of  the  helicopter  moiionr;  and  is  of  primary  importance  in 
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helicopter  stability  and  control.  The  term  "qua  si  "Static"  therefore  Implies 
Inclusion  of  pitch  rate  and  roll  rate  as  steady  state  rotor  variables,  and  in¬ 
stantaneous  adjustment  of  flapping  to  changes  of  those  variables. 

The  virtue  of  this  assumption  Is  the  reduction  of  the  flapping  differ¬ 
ential  equation  to  an  algebraic  equation  which  permits  elimination  of  flapping 
as  a  variable  in  the  equations  as  motion,  reducing  the  number  of  degrees  of 
freedom  and  the  characteristic  order  of  the  system. 

The  validity  of  this  assumption  is  primarily  dependent  on  the  fre¬ 
quency  of  the  helicopter  motion  n  ,  the  rotational  speed  of  the  rotor  &  , 

*50  n 

and  the  blade  Lock  number  y  ♦  It  is  snown  in  Reference  9  that  if 

y>ir 

then  this  assumption  introduces  little  or  no  error.  This  is  usually  the  case 
for  long  period  helicopter  motion. 

Although  the  assumption  as  described  has  been  widely  accepted  for 
Stability  and  Control  studies,  it  should  be  noted  that  the  dynamic  lag  of  cer¬ 
tain  stabilizing  devices  (the  analogue  of  the  bla^e  dynamic  lag  neglected  here) 
is  fundamental  to  their  operation,  and  cannot  be  disregarded  as  above.  It 
should  hardly  be  necessary  to  note  that  for  blade  flutter  or  strength  calcula¬ 
tions,  no  such  assumptions  are  permissible. 


3.4.  Rotor  Angular  Velocity 

The  rotor  angular  velocity  is  assumed  to  be  constant.  It  appears  that 
the  changes  in  rotor  torque  due  to  small  displacements  of  the  helicopter  from 
trim  are  typically  small  enough  that  the  resulting  rpm  changes  have  negligible 
effect  on  the  rotor  forces  and  moments.  This  assumption  is  fairly  well  justi¬ 
fied  by  test  data  (References  10  and  11)  in  which  rpm  changes  during  heli¬ 
copter  dynamic  response  are  negligibly  small.  The  simplification  obviates 
the  need  for  consideration  of  the  difficult  torque  equation,  and  the  evaluation 
of  the  torque  derivatives.  There  is,  however,  some  theoretical  evidence 
(Reference  12)  that  in  autorotation,  the  effects  may  be  more  pronounced.  The 
assumption  in  that  case  appears  conservative,  and  hence  may  ewn  there  be 
permissible. 
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3.5.  Coupling  of  Lateral  and  Longitudinal  Motions 

lateral  and  longitudinal  helicopter  motions  are  assumed  to  be  Inde¬ 
pendent — that  is,  lateral  motions  do  not  induce  or  affect  longitudinal  forces 
or  moments,  and  vice  versa.  This  assumption  is  one  which  definitely  limits 
the  following  theory  to  small  perturbations  of  aircraft  motion,  if  quantitative 
accuracy  is  expected.  Even  f>r  vanishingly  small  perturbations,  there  are 
some  such  coupling  effects,  but  these  are  generally  weak,  and  it  has  been 
found  that  reasonable  and  useful  results  can  be  achieved  without  their  con¬ 
sideration.  At  the  present  time  the  accuracy  of  estimation  of  major  stability 
derivatives  is  so  limited  that  consideration  of  relatively  weak  coupling,  with 
tremendous  complication  of  the  analysis,  is  not  justifiable.  As  the  small 
perturbation  theory  becomes  better  established,  interest  will  no  doubt  focus 
on  large  displacement  maneuvering  and  stability  problems,  and  coupling  will 
receive  further  study.  For  the  present  work,  the  assumption  obviates  the 
need  for  evaluation  of  many  so-called  "cross,  -  or  "coupling,"  derivatives. 

3.6.  Higher  Order  Terms 

In  the  formulas  quoted  in  the  following  sections,  an  attempt  has 
been  made  to  compromise  between  theoretical  accuracy  and  brevity,  parti¬ 
cularly  in  the  higher  order  fi  terms .  The  blade  stall  and  reversed  flow 
assumptions  severely  limit  the  theoretical  validity  at  large  advance  ratios, 
so  that  no  pretense  for  accuracy  is  made  above  about  /Is. 4  .In  view  of 
this,  higher  order  fl  terms  representing  less  than  about  20  percent  of  over¬ 
all  term  at  fls. 4  ,  are  omitted  from  the  formulae  entirely. 

It  should  be  emphasized  that  on  this  account,  and  also  because  of 
the  limiting  basic  assumptions ,  the  formulations  herein  cannot  be  relied  upon 
for  accuracy  (say  5  to  10  percent)  above  about  fl  =.3  ;  and  that  the  accuracy 
should  be  expected  to  deteriorate  rapidly  with  increase  in.  ft  ,  so  that  above 
fl  s.5  ,  even  trends  may  be  predicted  incorrectly. 
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4.  Rotor  Aerodynamic  Forces  and  Derivatives 
4.1.  Rotor  Thrust.  and  ^8  Forces 

The  most  general  statement  of  the  rotor  aerodynamic  forces  is  in 
terms  of  components  in  the  wind-  ANF  coordinate  system.  These  are  thrust, 
,  along  the  ANF  ;  Hw  force,  perpendicular  to  ANF  ,  in  the  plane 
of  the  wind  and  ANF  ;  and  force,  perpendicular  to  the  plane  of  the  wind 
and  ANF  .  The  forces  can  be  expressed  in  terms  of  the  basic  aerodynamic 
variables  measured  in  the  same  axis  system:  tJ  ,  the  advance  ratio;  9 , 

L  3  collective  pitch;  ,  the  inflow  ratio,  involving  ,  the  angle-of- 

attack,  measured  between  the  remote  wind  vector  and  a  plane  perpendicular 
to  ANF  ;  pitch  ard  roll  rates,  q  and  p  ;  and  blade  flapping  resolved  in 
the  same  axis  system,  given  by  Oq  ,  0|^  ,  and  b|^  . 

The  basic  equations  are 
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All  the  variables  in  the  above  are  taken  in  the  wind-  ANF  system. 
The  ones  for  which  the  distinction  is  necessary  are  ,  Oj^  ,  and  bj^  , 
which  are  so  identified  by  the  subscript.  The  equations  are  general,  since 
they  apply  whether  or  not  the  rotor  is  sideslipping,  pitching  or  rolling  about 
its  hub;  for  non-articulated  (rigid)  as  well  as  flapping  rotors;  regardless  of 
pitch-cone  or  pitch-flap  coupling  (as  with  8^  );  etc. 

For  articulated  rotors,  the  flapping  constants  in  the  above  are  given 
with  similar  generality  as 


I 
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The  average  inflow,  Xyy  ,  satisfies  the  equation 


•  n  ««* 


P) 


and  X  |  ,  the  fore-and-aft  induced  velocity  variation  is  given  by  Figure  2, 

X  yy  ,  Figure  1.  These  induced  flow  relations  aie  based  on  Reference  2,  and 
are  consistent  with  the  assumptions  concerning  induced  velocity  discussed 
earlier. 


The  above  equations  suffice  for  calculation  of  the  rotor  forces,  giver 

the  blade  characteristics  8  ,  0  ,  <T  ,  y  ,  Mg  ,  and  Ii  ;  and  the 

a  P  Q 

rotor  operating  conditions  U  ,  aw  ,9  ,  ~  ,  and  .  For  example, 

I  M  U 

one  might  solve  for  Cy  ana  '  w  by  trial  and  error  using  (1)  and  (7);  then 

W 

Oq  ,  0|w  ,  and  bjw  ;  ?  d  finally,  Cj^  and  Cy^  .  This  procedure 

would  not  usually  be  convenient,  however,  for  several  reasons: 

1)  the  operating  variables,  9  and  CKyy  ,  are  not  usually  known  at 
the  beginning,  but  must  rather  be  found  from  trim,  or  equilibrium,  conditions 
involving  Cy  ana  CH  directly; 

2)  with  pitch -flap  coupling,  as  with  8^  ,  the  ANF  position  varies 
with  flapping,  and  hence  dyy  is  not  known  originally  but  depends  on  CJ| 
and  b  |  ; 


3)  the  process  of  taking  derivatives  of  the  forces  with  respect  to  the 
operating  variables  is  unduly  complicated;  and 

4)  it  would  be  more  convenient  to  have  the  forces  resolved  along  and 
perpendicular  to  the  rotor  shaft,  in  directions  fixed  in  the  helicopter,  rather 
than  moving  with  control  deflection,  sideslip,  flapping,  etc. 

These  objections  can  largely  be  removed  by  eliminating  8  in  the  H 
and  Y  force  equations  in  favor  of  T  ,  by  using  equation  1;  by  substituting 
equations  4,5,  and  6,  for  flapping,  into  equations  1,  2,  and  3  for  the  forces; 
and  by  forming  new  expressions  for  the  forces  along  shaft  axes  in  terms  of 
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aircraft  operating  conditions .  The  last  is  accomplished  by  means  of  the  axis 
transformation  equations  given  below: 


T*  *  T* 

(8) 

Y,  »  Yw  -  HwSl  +  Tw  A„ 

(9) 

**  “  Tw  B|# 

(10) 

where 

0i“0s  +  -Y'  +  Q 

(11) 

8  q 

®wcs  ®s”  ®ls  “  v 

(12) 

blw~  bls  “Als+  +Bls) 

(13) 

°lw=  °ls  +  Bls  ■  ^i(bl«  _Als) 

(14) 

In  the  above  equations  the  ANF  directions  relative  to  the  rotor 

shaft  appear  for  the  first  time .  They  are  A|s  and  , 

and  depend  on 

cyclic  pitch  control  settings  (through  the  swashplate)  and  on  flapping  with 

respect  to  the  shaft,  if  pitch -flap  coupling  exists,  as  with 

8a  .  With 

pitch-flap  coupling  given  by  the  symbol  T  ,  the  blade  pitch  and  flapping 

are  related  by: 
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The  proceeding  sets  of  equations  suffice  to  find,  by  successive  elimina¬ 
tion  of  unwanted  variables,  expressions  for  the  rotor  forces  and  their  derivatives, 
resolved  along  shaft  axes,  in  terms  of  basic  rotor  operating  conditions.  These 
expressions,  along  with  those  for  the  derivatives,  are  tabulated  in  tables  1 
through  3.  They  are  given  for  equilibrium,  steady,  symmetrical  flight  condi¬ 
tions  ,  for  which  (p)j  *  (q)j  ■(/aji-(+)i’0. 

Formulas  in  the  left-hand  column  of  the  tables  are  general,  in  that  they 
apply  for  cases  involving  pitch-flap  coupling,  and  can  be  used  for  any  (7  up 
to  roughly  jJ  a  .*5  .  They  involve  a  series  of  "  f  "  factors  which  depend  on 

Lock's  constant,  y  ,  and  the  pitch-flap  coupling  parameter,  T  .  The  *  f  " 
factors  can  conveniently  be  found  from  the  charts,  Figures  8  to  19.  For  the 
case  of  no  pitch-flap  coupling,  these  "  +  "  factors  reduce  to  simple  constants, 
many  of  them  equal  to  zero.  For  the  Hohenemser  "pitch-cone"  rotor,  which 
couples  collective  pitch  and  coning,  but  not  cyclic  pitch  with  cyclic  flapping, 

T  should  be  taken  zero  for  all  the  "  f  "  factors  except  fjj  ,  fg  ,  fy  , 
and  “fg 

One  of  the  rotor  variables  is  the  wind  axis  inflow  ratio,  X^  ,  defined 
by  equation  7.  It  has  not  been  eliminated  because  of  the  algebraic  difficulities 
of  this  equation.  The  chart,  Figure  1,  has  however,  been  prepared  as  a  con¬ 
venient  means  to  find  Xw  as  a  function  of  ,  fl  ,  and  Cy  ,  for  use  in 
the  force  equations.  For  the  same  reasons,  the  force  derivatives  are  in  terms 
of  partial  derivatives  of  X^  .  These  have  been  calculated  from  equation  7 
and  can  easily  be  evaluated  from  the  charts,  Figures  3  to  4. 
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In  the  case  of  pitch -flap  coupling,  the  ANF  angle-of-attack ,  a  ^  , 
of  which  X^is  a  function,  can  be  given  as 

a*--!TT?-c+8a{(4?),(t;  -  f2 

“i|Xw(+3  -+4°2)}  (l8) 


where 


t„CT)-vo) 


The  above  relation  is  simpler  for  the  pitch-cone  case  and  the  case  of 


no  pitch -flap  coupling  at  all,  since  it  reduces  to 


«w  *  “s  -  B|c 


The  general  formulas  of  the  tables  can  sometimes  be  further  simplified, 
as  follows: 

1)  For  fi  less  than  about  .15,  the  terms  can  safely  be  neglected. 
A  very  little  experience  with  the  formulas  will  indicate  which  terms  are  negligi¬ 
ble  for  the  type  rotor  being  considered. 

2)  For  u  greater  than  about  .15,  the  charts  for  Xu  and  its  partial 

2  2 

derivatives  are  unnecessary,  since  it  can  be  assumed  that  X  «  fi 
and  therefore 


.  .  Cj 

*w  ®w  "  2p. 
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Cj 

°w  +  2^“ 


The  maximum  possible  simplifications  of  the  general  formulae  are  for  the  no 
pitch-flap  coupling  case  (r.O)  at  advance  ratios  ver  near  .15.  At  this 
fi  ,  the  higher  order  fi  terms  can  be  dropped,  and  eliminated  by  the 
above  approximations,  yielding  the  "special  case"  formulas  in  the  right-hand 
columns  of  the  tables.  These  are  included  for  their  simplicity  and  value  for 
preliminary  estimates,  but  the  restrictions  on  their  validity  should  not  be  over¬ 
looked. 

Several  of  the  general  formulas  for  the  rotor  forces  and  derivatives 
contain,  for  brevity,  lumped  parameters,  the  formulas  for  which  are  collected 
in  Table  6,  for  easy  reference. 

4.2.  Trim  Calculate  ;s 

The  expressions  for  the  rotor  forces  and  moments  and  their  deriva¬ 
tives  of  Tables  1  through  3  contain,  as  a  basic  parameter,  the  inflow  ratio, 
Xw.  It  has  not  been  eliminated,  as  previously  discussed,  because  of  the 
algebraic  complexity  of  equation  7.  Before  the  forces  and  derivatives  can  be 
calculated  for  any  flight  condition,  the  must  be  determined.  The  prob¬ 
lem  is  to  find  a  value  for  Xyy  for  which  the  forces  acting  on  the  helicopter 
are  in  equilibrium  for  the  flight  condition  specified. 

The  process  inevitably  involves  consideration  of  lift,  drag,  and 
pitching  moment  equilibrium  of  the  aircraft,  and  can  be  quite  an  involved 
iteration  procedure  unless  advantage  is  taken  of  charts  prepared  specifically 
for  the  purpose.  The  problem  is  so  basic,  however,  that  it  seems  appropriate 
to  discuss  at  least  one  possible  technique  at  this  point. 
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One  can  usually  at  the  outset  calculate,  for  the  specified  flight  condition, 
the  rotor  thrust  coefficient  required  for  lift  equilibrium;  and  also  the  pitching 
moment  of  the  rotor  forces  about  the  aircraft  center  of  gravity  required  to  balance 
the  fuselage  and  tail  moments.  The  latter  condition  will  specify  the  direction 
of  the  rotor  force  vector,  Qg  .  The  procedure  of  satisfying  equilibrium  of  drag 
forces  will  depend  on  whether  the  flight  condition  specifies  climb  angle  or 
engine  power.  Either  would  determine  Xy^and  hence  the  other.  Both  cases 
will  be  discussed  below. 

A  convenient  parameter  is  h4  ,  defined  as  HQfor  7*0.  h4  is  a  function 
of  C  and  X^given  in  equation  79  and  displayed  in  Figure  5.  Using  equa¬ 
tions  18  and  45,  it  may  be  shown  that  »  19  r  \ 

u 

"4  ZCj  (19) 

a  <t 

Now  drag  equilibrium  of  the  helicopter  requires  that 

,  Df 

°s  +  as+  ^  *w“  (20) 

Consider  first  the  case  where  climb  angle,  «  is  specified.  From  (20), 
(dg+flg)may  be  found.  Then  to  find  <tyy,  assume,  say,  three  values  for  it  and 
find  corresponding  va’ues  of  H4  from  (19),  and  X  from  Figure  1.  Plot  these 
pairs  of  H4'X*  on  Figure  5.  The  correct  value  of  X^ is  the  intersection  of 
this  line  with  the  appropriate  fl  line  of  the  lower  family  on  the  chart.  The 
rotor  forces,  their  derivatives,  and  torque  can  now  be  determined.  Since  dg 
is  known  from  moment  equilibrium,  dgis  known,  and  the  longitudinal  control 
to  trim  follows  from  eqi  tion  18.  If  necessary,  the  assumed  average  blade 
drag  coefficient,  8  ,  may  be  improved  by  using  the  charts  for  the  torque 
equation,  as  discussed  in  the  next  section. 

For  X)  greater  than  about  .15,  the  above  procedure  can  be  greatly 
simplified  by  using  the  parameter 

h;-h4+^ 


shown  in  the  upper  part  of  Figure  5.  It  can  be  shown,  using  the  approximate 
expression  for  Xyywith  (19),  that 
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in  which  <X^  does  not  appear.  For  cases  of  U  >.|5  *  may  be  calc u 
lated  by  (21),  and  X^  follows  immediately  from  the  upper  set  in  Figure  5. 

X  then  is  known,  and  dyy  may  be  found  from  Figure  1  or  equation  7.  The 
other  unknowns  can  thereafter  be  calculated  in  the  same  way  previously 
described. 

In  the  other  event,  where  rotor  torque  is  specified  by  the  flight  condi¬ 
tion  but  the  climb  angle  is  unknown,  X  is  found  first  from  charts  for  the 
torque  equation,  as  given  in  the  next  section.  At  the  same  time,  tho  mean 
effective  blade  drag  coefficient,  8  is  evaluated.  Knowing  X  ,  Uy  may 
be  found  from  Figure  1,  and  then(o'g  4* (^follows  from  or  using 
Figure  5  and  (19)  or  (21).  Finally,  from  (20),  the  climb  angle  «  may  be 
calculated. 

4.3.  Rotor  Torque 

If  the  engine  power  absorbed  by  the  main  rotor  is  known,  then  a  dif¬ 
ferent  method  is  used  to  determine  the  inflow  ratio  Xyy  .  Charts  are  pre¬ 
sented  for  the  torque  equilibrium  of  a  rotor  from  equations  given  in  Refer¬ 
ence  13,  with  the  collective  pitch  angle  eliminated  in  favor  of  the  thrust 
coefficient,  and  the  expression  for  accelerating  torque  modified  to  account 
for  the  reverse  flow  region  in  the  manner  discussed  in  Section  3.2. ,  assump¬ 
tion  iv. 

The  equation  for  the  torque  coefficient  of  a  rotor  may  then  be  written 
as: 


(22) 
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where  Qj  «  #  Og  3  t(/*,X/^)and  Q^,5  f  (/*) 

The  complete  expression  is  given  in  table  6.  Normally  the  term  Qij  may  be 
neglected  as  small.  If  is  neglected,  then  A  may  he  determined  in  the 
fo llowing  way .  _ 

ZCt  a 

A*  the  outset,  ■  *  will  normally  be  known,  and  Qi  may  be  ob- 

Q  <7  I 

tained  from  Figure  6.  02  may  then  be  found  from  equation  22  (putting 

Qj»  0  )  and  thereafter  X ^  is  determined  by  Figure  7.  The  remaining 
trim  quantities  are  then  found  from  Figure  5  and  the  previous  procedure.  Of 
course,  the  process  may  be  reversed,  so  that  if  trim  conditions  are  deter¬ 
mined  for  a  given  climb  angle,  then  Figures  6  and  7  may  be  used  to  deter¬ 
mine  the  torque  coefficient. 

The  small  term  involving  Qj  can  easily  be  included  by  a  simple 
iteration.  After  initially  neglecting  03  to  find  X  ,  then  compute  the 
term  and  include  it  in  equation  22  to  correct  q2  .  A  corrected  X  ^  follows 
from  Figure  7.  The  term  is  so  small  that  a  single  iteration  surely  will 
suffice. 

An  average  blade  drag  coefficient  for  use  in  the  rotor  force  equations, 
may  be  determined  from  the  equation 


where  Cq  is  determined  as  described  above,  and  the  quantity  enclosed  in 
brackets  may  be  found  directly  from  Figure  5  for  t  .15  ,  for  fl  <.!5  , 
h4  may  be  determined  from  the  same  figure  and  the  additional  term  com¬ 
puted  directly. 
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5.  Rotor  Hub  Moments 

In  addition  to  the  aerodynamic  forces  developed  by  a  rotor,  the 
aerodynamic  and  inertial  forces  acting  on  an  articulated  rotor  blade  pro¬ 
duce  moments  about  the  rotor  hub  when  the  flapping  hinges  are  not  located 
on  the  rotor  shaft.  These  moments  arise  from  the  vertical  shear  force  on 
the  flapping  pin  acting  at  a  distance  t  ,  the  eccentricity  of  the  flapping 
hinge,  from  the  rotor  shaft.  The  constant  part  of  the  shear  force  is  the 
rotor  thrust  and  does  not  produce  moments.  However,  the  cosine  and  sine 
components  produce  pitching  and  rolling  moments,  respectively,  about  the 
rotor  hub.  In  the  wind-  ANF  system  the  rotor  hub  moments  may  be  ex¬ 
pressed  as: 
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As  discussed  in  section  3.4.,  it  is  more-convenient  to  resolve  rotor 
forces  and  moments  in  the  shaft  axis  system.  The  hub  moments  in  this  axis 
system  may  be  expressed  in  terms  ■'ji  the  wind**  ANF  quantities  as: 

M$  -  Mw  +  Lw  $t  (24) 

Ls  *  Lw  «  Mw  0S  (25) 

In  tables  4  and  5  these  quantities  and  their  derivatives  are  presented 
in  terms  of  the  basic  rotor  variables  Cy  #  X  ,  and  fi  .  These  tables  are 
in  the  same  form  as  the  rotor  force  tables,  that  is,  the  general  equation  is 
presented  on  the  left-hand  side,  and  the  simplified  form  for  fi  .15  is 
presented  on  the  right.  The  discussion  of  section  3.4.  concerning  the  appli¬ 
cability  of  these  formulae  holds  true  here. 
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TABLE  I  -  ROTCR  THRUST  AND  ITS  DERIVATIVES 


TABLE  2 -ROTOR  H  FORCE  AND  ITS  DERIVATIVES 
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TABLE  4  ~  ROTOR  HUB  PITCHING  MOMENT  AND  ITS  DERIVATIVES 


TABLE  5 -ROTOR  HUB  ROLLING  MOMENT  AND  ITS  DERIVATIVES 


EQU. 

NO. 

GENERAL  FORMULA 

SPECIAL  CASE  FORMULA  FOR  : 
r«  0  .15 
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dCf  j  » 
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4cr,  dx 
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TABLE  6-  INFLOW  DERIVATIVES  AND  AUXILIARY  PARAMETERS 
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FIG.  I  INFLOW  RATIO  VS  ADVANCE  RATIO 
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FIG.  4  INFLOW  RATIO  PARTIAL  DERIVATIVES 
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FIG.  7  TOR^JE  EQUATION  FACTOR  ,  Q2 
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FIG.  9  THRUST  EQUATION  FACTORS,  5,  6,  7,  8 
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FIG.  13  Y  FORCE  EQUATION  FACTORS,  17,  18 
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FIG.  17  HUB  MOMENT  EQUATION  FACTORS ,  26,27 


83  PARAMETER,  r 

FIG.  18  HUB  MOMENT  EQUATION  FACTORS,  28,29,30 
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FIG.  19  THRUST  DERIVATIVE  FACTORS,  31,32 


